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Abstract 

In the framework of the minimal SU(3)c ® SU(3)l ® U(1)jv model, the lepton-flavour- 
violating decay ty + — > fi~ e + e + is calculated without directly invoking lepton mixing. The 
branching ratio for this rare pion decay mode is found to be much smaller than the current 
' experimental upper limit. Dropping out anomalous interactions, this result coincides with 

the previous calculation . 
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At present, neutrinos are presumably massive and mixed as indicated in various exper- 
^| iments: SuperKamiokande |IJ and others 0. This significant deviation from the standard 

■ model (SM) calls for its extension. The models based on the SU(3)c ® SU(3)i (g) U(l)jy (3 

3 1) gauge group |3|,|4| are one of the most popular in such extensions beyond SM. The SM 
assumes lepton-flavour-number conservation, and its observed violation would be a clear in- 
dication of new physics. In the 3 3 1 models the lepton-flavour number is not conserved, and 
these models have motivated a variety of dedicated sensitive searches for rare decay modes 
of muons and kaons and for neutrino oscillations ||. It is known that the muon system is 
one of the best places to search for lepton flavour violation, compared with the others. The 
"wrong" muon decay fjT — > e'Vgp^ is widely used to put a lower bound on the singly-charged 
bilepton mass (M Y > 230 GeV) §. 

In this work we pay attention to the lepton-flavour- violating pion decay tt + — > /i~ e + e + . 
The upper limit in its branching ratio is given R < 1.6 x 10~ 6 at 90 % confidence level 
By suggesting the lepton mixing or horizontal interactions, the above decay has been studied 
theoretically in Ref. M. However, this decay may be described by the minimal 3 3 1 model 
in simple manner without directly invoking lepton mixing. 

To start, we firstly give some basic elements of the model (for more details see [0). Three 
lepton components of each family are in one triplet, 
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where a = 1, 2, 3 is the family index. Under SU(3)l, two of the three quark families 
transform as antitriplet and one family transforms as triplet, 





QiL = -u iL , {i = 1,2), Q 3L 



The right-handed quarks are singlets under SU(3)l- The exotic quarks T and D{ have an 
electric charge + 5/3 and - 4/3, respectively 

There are five new gauge bosons: the Z' and the charged bileptons with lepton number 
L = ±2, which are identified as follows: Y~ = W* - iWj, y/2 X~ = Wj - zWj, and 
their couplings to leptons are given by |TT| : 

qCC _ 9 



2V2 



*f{i - i,)cFy- - h^cFx-- + h.c. . (1) 



The interactions among the charged vector fields with quarks are 

+ (u 3L YT L - D iLl »d iL )Y- + h.c.]. (2) 

It is to be noted that the vector currents coupled to X , X ++ vanish due to Fermi statistics, 
and the exotic quarks interact with ordinary ones only via the bileptons and non-SM Higgs 
bosons. 



The current experimental lower bound on the exotic quark mass is 200 GeV ||12j1 , while 
the lower bound on the bilepton mass is in the range of 300 GeV. 

To deal with the above process we also need the coupling constants of the bileptons X, Y 
to the SM weak- vector boson W. In the notation of Refs. [13| it is: CWXY = 

Now we start with the decay 

n +( K ) - /i-(p) u„{q) e + (h) e + (k 2 ), (3) 

where the letters in parentheses stand for the momenta of the particles. We assume that the 
Higgs bosons responsible for lepton-flavour- violating interactions as well as the exotic quarks 
are much heavier than the standard model W boson. Hence the contributions from the exotic 
quarks and non SM Higgs bosons are suppressed. With new gauge bosons carrying lepton- 
number L = 2, the process (Q) can be described simply by the Feynman diagram depicted 
in Fig. 1. 

For small momentum-transfer (q 2 << rriyy , M\ , My) , as is the case here, the matrix 
element for this process is found to be 



Mfi = 2^i^[-(P + K) p K, + (K + L),K,+ (-L + P).Kg, y ] 

xu^h^l - ^Cu^.vJih)^^^), (4) 
where the following combinations of four vectors are introduced 

P = h + k 2 , Q = h - k 2 , L = p + q, N = p - q, K = P + L. (5) 
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Figure 1: Feynman diagram for the decay n + (K) — > /i (p) ^ t (g) e + (A; 1 ) e + (A;2) 

in the 3 3 1 model 



The squared matrix element is given by 



\M 



fi\ 



128 



x 

X 



My My 



PmQn 



kjkt + kt kl - grt \k x .k 2 - m 



(6) 



where the notation C/3 7 = [— (P + K)pK 7 + (K + L) 7 Kp + (— L + P).iT g^] is used. 

In order to describe the kinematics of the decay, we introduce the following vectors: Let v 
be a unit vector along the direction of flight of the dipositron in the ir + rest system (£„-), a 
be a unit vector along the projection of the three-momentum of the e + in the e + e + center- 
of-mass system (S 2e ) perpendicular to and b be a unit vector along the projection of the 
three-momentum of the in the fj,~ center-of-mass system (E^) perpendicular to —v. 
Then the kinematics of this decay is similar to the one given in jTJ], which consists of five 
variables: s e = P 2 , = L 2 , and three angles: (i) e , the angle of the e + in E 2e with respect 
to the dipositron line of flight in £„-, (ii) 0^, the angle of the [i~ in T,^ with respect to the 
[j,v n line of flight in H n , and (iii) (ft, the angle between the plane formed by the positrons in 
and the corresponding plane formed by the //~, v^. The angles e , 6^ and are shown in 
Fig. 2. 




Figure 2: Illustration of the angles e , 0^ and 
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Then the decay width for the pion decay (|3]) is written as 



dT = ( o ) 914 6 3 51 l-^/il 2 (1 - z n) cr e X ds e ds /1 d(cos e )d(cos fl )d(f). (7) 



7r spins 



In (0), (|J is the statistical factor indicating that two (identical) positrons in the final 
state ||15|| . With the above definitions we have the following scalar products 

Q 2 = 4m 2 - Se , N 2 = 2ml ~ s ^ K * = m ^ L - N = m l> 

P.L = h m 2 e -8 e -8,i), P.N = z u P.L+ {I- Zfl )X cos 9 U , Q.L = a e X cos9 e , (8) 

Q.N = z^Q.L + er e (l — z^)P.L cos 9 e cos 9^ — (s^) 1 / 2 sin 9 e sin 9^ cos <f>, 
d = E^L u N„P a Q(3 = -(s e s^ 2 a e (l - z fl )Xsm0 e sm0 tl sm<f>, 



where 



z, = ^, a e ee ( 1 - ^V /2 , X ee ((P.L) 2 - s e s,y/ 2 , 



and m e , m M , m„ stand for masses of the electron, the muon and the pion, respectively. 
The range of the variables is 

Am 2 < s e < (uin — m Ai 



■m 



2 / „ / /^\2 



< -V < (m w - , (9) 
< e , Op < 7T, < (f> < 2tt. 

It is to be noted that an imaginary part of |A4/j| 2 connected with pseudotensor d is linear 
in sin</>, i.e. no such a term like Q.Nd, hence it will be removed after integration over the 
angle <p. In resulting we get the decay width being a real number, as it has to be. 

The integrations over the angles can be carried out analytically by using Mathematica. 
The numerical integrations over the effective masses squared s e and s M are carried out by 
employing the Monte Carlo routine VEGAS [16]. After changing to dimensionless parameters 
x e = Vn = we get the decay width 

where N is numerically evaluated, iV = 6.17 x 1CT 6 . We recall that the main (99.987 %) 
decay mode of the n + is well-known 

n 2 f 2 m 2 

T(n + -> v,) = 5 3 (m* ~ ml) 2 ~ 2.63 x 1(T 17 GeV. (11) 



From and ([TT]) we get the branching ratio 



R-n 



r(7r+ -> z/ M e+ e+) 6.17 x 10" 6 G 2 F m^m 



i-i 



r(vr+ -> ^) 327r 5 M£ M*ml(ml - ml) 2 

U7xir "«ppp^' (12 » 



Putting Mx ~ My ~ 120 GeV as a lower limit obtained from the LEP data analysis ||17|| , 
we get R-x ~ 2.3 x 1CT 34 . This number is much smaller than the current experimental 
upper limit, but it coincides with the previous theoretical evaluation without anomalous 
interactions included |§. It rises a question about the mechanism for large lepton-flavour- 
violating pion decay mode. However, it is worth mentioning that the experimental data on 
R n decrease with time, for example the 1988 data were < 8 x 1CT 6 , while the 1998 data 
are R n < 1.6 x 1CT 6 . We suggest that by adding contributions from diagrams with the exotic 
quarks and Higgs bosons the situation will be modified but not improved too much. 

Our calculation can be analogously applied for the lepton-flavour-violating kaon decay 
K + — > fi~h>^e + e + , which has an experimental branching ratio of Rk < 2.0 x 10~ 8 . However, 
the main decay mode K + — > /i + has only a branching ratio of 69.51 %, instead of 99.987% 
in the 7r + case considered here . 

In summary, we have considered the lepton-flavour-violating pion decay without directly 
invoking lepton mixing. Our result is by twenty eight orders smaller than the current ex- 
perimental upper limit. This conclusion should not be modified too much by including 
contributions from the exotic quarks and Higgs bosons. Hence, the mechanism for large 
lepton-flavour-violating pion decay mode is a mystery. 
Acknowledgments 

This work was supported by the DAAD grant. The author thanks Prof. D. Schildknecht 
for stimulating remarks and S. Dittmaier for consultation on the Monte Carlo routine VE- 
GAS. This work was initiated when the author was at LAPTH, Annecy, France. He thanks 
G. Belanger and F. Boudjema for consultation on the decay kinematics and a numerical 
calculation. Thanks are due to P. Aurenche and Theory Group at LAPP and Bielefeld 
University for kind hospitality. 



References 

[1] Y. Fukuda et ai, Phys. Lett. B433, 9 (1998); Phys. Rev. Lett. 81, 1562 (1998); T. 
Kajita, in Proceedings of the XVIIIth International Conference on Neutrino Physics 
and Astrophysics, Takayama, Japan (June 1998). 



[2] For the review see, S. M. Bilenky, C. Giunti, and W. Grimus, [|hep-ph/981236qi , Prog. 
Part. Nucl. Phys. 43, 1 (1999). 

[3] F. Pisano and V. Pleitez, Phys. Rev. D 46, 410 (1992); P. H. Frampton, Phys. Rev. 
Lett. 69, 2889 (1992); R. Foot, O.F. Hernandez, F. Pisano, and V. Pleitez, Phys. Rev. 
D 47, 4158 (1993). 



5 



[4] R. Foot, H. N. Long, and Tuan A. Tran, Phys. Rev. D 50, R34 (1994); H. N. Long, ibid. 
D 54, 4691 (1996); J. C. Montero, F. Pisano, and V. Pleitez, ibid. D 47, 2918 (1993). 

[5] M. M. Guzzo et al, IPT-P.057/99, ||hep-ph/99U8^DE| . 

[6] E. D. Carlson and P. Frampton, Phys. Lett. B283, 123 (1992); H. Fujii, S. Naka- 



mura, and K. Sasaki, Phys. Lett. B299, 342 (1993); H. N. Long and T. Inami, [|hep- 



ph/ 9902475| , to appear in Phys. Rev. D; see also the review Y. Kuno and Y. Okada, 



KEK-TH 639, [|iep-ph/9909265| . 

[7] Particle Data Group, R. M. Barnett, et al, Eur. Phys. J. C3, 1 (1998); Phys. Lett. 
B204, (1988). 

[8] V. A. Baranov, et al, Sov. J. Nucl. Phys. 54, 790 (1991); Yad. Fiz. 54, 1298 (1991). 

[9] V. D. Laptev, in Program of Experimental Inverstigations at the Meson Factory of the 
Institute of Nuclear Research of the USSR Academy of Sciences [in Russian] (Moscow, 
1984), p. 145. 

[10] D. Ng, Phys. Rev. D 49, 4805 (1994). 

[11] P. H. Frampton and D. Ng, Phys. Rev. D 45, 4240 (1992). 

[12] P. Das, P. Jain, and D. W. Mckay, |hep-ph/9808"25^1 , Phys. Rev. D 59, 055011 (1999). 

[13] T. Ishikawa, T. Kaneko, K. Kato, S. Kawabata, Y. Shimizu, and H. Tanaka, GRACE 
manual, KEK Report 92-19 (1993); H. N. Long and T. Inami, in [§. 

[14] A. Pais and S. B.Treiman, Phys. Rev. 168, 1858 (1968); N. Cabibbo and A. Maksy- 
movicz, ibid. 137, B438 (1965); erratum ibid 168, 1926 (1968); J. Bijnens, Nucl. Phys. 
B337, 635 (1990). 

[15] C. Quigg, Gauge theories of the strong, weak and electromagnetic interactions, Benjamin 
(1983); T-P. Cheng and L-F. Li, Gauge theory of elementary particle physics, Oxfort 
University press, New York (1984); Q. Ho-Kim and X. Y. Pham, Elementary particles 
and their interactions. Concepts and phenomena, Springer, Berlin Heidelberg (1998). 

[16] CP. Lepage, J. Comp. Phys. 27, 192 (1978). 

[17] P. Frampton, Int. J. Mod. Phys. A13, 2345 (1998). 



6 



